Domain Growth in Ising Systems with Quenched Disorder 
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We present results from extensive Monte Carlo (MC) simulations of domain growth in ferromag- 
nets and binary mixtures with quenched disorder. These are modeled by the random-bond Ising 
model and the dilute Ising model with either nonconserved (Glauber) spin-flip kinetics or conserved 
(Kawasaki) spin-exchange kinetics. In all cases, our MC results are consistent with power-law growth 
with an exponent 9(T, e) which depends on the quench temperature T and the disorder amplitude 
e. Such exponents arise naturally when the coarsening domains are trapped by energy barriers 
which grow logarithmically with the domain size. Our MC results show excellent agreement with 
the predicted dependence of 0(T, e). 
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I. INTRODUCTION 



Consider a binary mixture which is homogeneous at 
high temperatures. This system becomes thermodynam- 
ically unstable if it is quenched below the critical temper- 
ature. The subsequent evolution of the system is charac- 
terized by the formation and growth of domains enriched 
in either component. These domains have a characteris- 
tic size R(t), which grows with time. The domain growth 
law [R(t) vs. t] depends on general system properties, 
e.g., the nature of conservation laws governing the or- 
der parameter evolution; the presence of hydrodynamic 
velocity fields; the presence of quenched or annealed dis- 
order, etc. There is a good understanding of the growth 
laws for pure and isotropic systems 0, Q • For the 
case with nonconserved order parameter, e.g., ordering of 
a magnet into up and down phases, the system obeys the 
Lifshitz-Cahn- Allen (LCA) growth law, R(t) ~ i 1/2 . For 
the case with conserved order parameter, e.g., diffusion- 
driven phase separation of an AB mixture into A-rich 
and B-rich phases, the system obeys the Lifshitz-Slyozov 
(LS) growth law, R(t) ~ t 1 / 3 . 

Recent interest in domain growth problems has focused 
on modeling and understanding the effects of various ex- 
perimentally relevant features. In this context, an im- 
portant set of analytical and numerical studies has in- 
vestigated coarsening i n sy stems with quenched disorder 

In general, one expects 
that trapping of domain boundaries by disorder sites will 
result in slower domain growth. However, these studies 
were unable to clarify the nature (or even existence) of a 
universal growth law. In a recent letter [T^ , we have re- 
visited this problem through comprehensive Monte Carlo 
(MC) simulations of kinetic Ising models. In our letter, 
we presented MC results for ordering in random mag- 
nets, modeled by the random-bond Ising model (RBIM) 
with nonconserved (Glauber) spin-flip kinetics. (In the 
RBIM, the presence of disorder is mimicked by randomiz- 
ing the exchange coupling between spins.) In this paper, 
we present further results for coarsening in two classes of 
disordered systems: 



(a) The RBIM with conserved (Kawasaki) spin-exchange 
kinetics, which models phase separation in disordered bi- 
nary mixtures. 

(b) The dilute Ising model (DIM) with both nonconserved 
and conserved kinetics. The DIM is relevant in cases 
where disorder is introduced via either bond dilution or 
site dilution. 

The results in this paper, in conjunction with those in 
our letter, constitute a novel understanding of domain 
growth in systems with quenched disorder. This paper is 
organized as follows. In Sec. II, we summarize arguments 
for growth laws in disordered systems. In Sec. Ill, we 
present results for the RBIM with conserved kinetics. In 
Sec. IV, we present results for the DIM with both noncon- 
served and conserved kinetics. Finally, Sec. V concludes 
this paper with a summary and discussion of our results. 



II. GROWTH LAWS IN DISORDERED 

SYSTEMS 

A. Nonconserved Case 

An important step towards understanding growth laws 
in nonconserved systems is due to Lai et al. (LMV) 0] . 
LMV proposed four classes of systems, determined by 
the dependence of the energy barrier to coarsening on 
the characteristic scale. The growth of domains is driven 
by a curvature-reduction mechanism as 



dR _ a(R, T) 
~dt ~ R : 



(1) 



where the diffusion constant a(R,T) depends on the 
domain scale R and temperature T, in general. For 
pure systems, the diffusion constant is independent of 
the length scale, i.e., a(R,T) = ao- The corresponding 
growth law is the LCA law, R(t) = {2a t) 1/2 . 

Let us next consider systems with quenched disorder. 
At early times and small length-scales, the growing do- 
mains are not affected by disorder [a(R,T) ~ ao] and 
the growth law is the same as that for the pure case. 
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At late times, the domains are trapped by disorder sites, 
creating a barrier (Eg) to domain growth. Then, the 
asymptotic dynamics is driven by thermal activation over 
disorder barriers with a(R,t) ~ ao exp(— (3Eb), where 
(3 = T _1 (ks = lj. For the bond-disordered case, Huse 
and Henley (HH) y| argued that the energy barrier scales 
as Eb(R) — eR^, where e is the disorder strength. The 
barrier exponent ip depends on the roughening exponent 
£ and the pinning exponent x as V 1 — x/(2 — 0- Fur- 
ther, the roughening and pinning exponents are related 
as \ = 2( + d — 3, where d is the dimensionality. For 
power-law barriers, Eq. yields an asymptotic growth 
law which is logarithmic, viz., 



R(t) 
to 



1 

aoip 



2/V 



(2) 



We can reformulate the early-time and late-time behav- 
iors as limiting cases of a crossover function: 



/?(/) = RtAT.oh [~ ) 



where 



and 



Ro(T, e) = 



T 



1/2 



(3) 



(4) 



(lnz) 1 ^. 



a; < 1, 
a: » 1. 



(5) 



For d = 2, C = 2/3 and X = 1/3 [H H3, yielding 
ip = 1/4. For e? = 3, a perturbative calculation gives 
ip ~ 0.55 |5j. There have been a number of numerical 
simulations]!, IS H OH IH Ed and experiments 
[TsL TTgL l20| which have attempted to test the HH scenario. 
However, to date, there is no clear confirmation of HH 
growth in the asymptotic regime. As a matter of fact, it 
is not even clear whether there is a universal law which 
characterizes the disorder-affected growth regime. 

In recent work |14|. we have revisited this problem 
via extensive MC simulations of the RBIM with non- 
conserved kinetics. Our results were consistent with 
power-law domain growth, but with a temperature- and 
disorder-dependent exponent. Similar observations have 
been made in experiments on coarsening in disordered 
systems |0, E3 • Such growth exponents can be un- 
derstood in the framework of a logarithmic (rather than 
power-law) independence of trapping barriers. In the 
context of the DIM, Henley [2l| and Rammal and Benoit 
[22| have argued that the fractal nature of domain bound- 
aries results in a logarithmic independence of energy bar- 
riers. We propose that this is generally applicable and 



examine the implications thereof [23j. Recall that, at 
early times and small length scales, we expect disorder- 
free domain growth. Then, the appropriate logarithmic 
barrier-scaling form is as follows: 



E B (R)^eln(l + R), 



(6) 



where R is measured in dimensionless units. Substituting 
a(R,T) ~ a exp(— /3£b) in Eq. (JTJ, we obtain 



dR 
~dt 



R 



(l + RY 



e/T 



(7) 



The solution of Eq. Q is 
R{t) ~ (2a t) 1/2 , 

( 2 + ?) 



t < to, 

1 e(T,e) 

dot 



t > to, 



with the asymptotic growth exponent 
6{T,e)= ' 



e/T 



(8) 



(9) 



The crossover length and time can be identified by rewrit- 
ing Eq. JSJl in the form of Eq. with 



and 



Ro 



to = 



h(x) 



1 



(20)9/(1-29)' 
1 1 

~(26> 2e ) 1 /(i-20)' 



(10) 



l '\ z«l, 
e , ir>l. 



(11) 



In our letter (l3|. we have shown that the growth ex- 
ponent for the nonconserved RBIM is consistent with 
Eq. ©. Let us next discuss the implications of power- 
law and logarithmic barriers for domain growth with con- 
served kinetics. 



B. Conserved Case 

In the absence of disorder, the domain scale obeys the 
Huse equation |24j : 



dR 
~dt 



Do 
R 2 '' 



(12) 



with the solution R(t) = (SDot) 1 ^ 3 . The presence of 
disorder renormalizes the diffusion constant Dq by an 
Arrhenius factor: D(R,T) ~ D exp(— (3Eb)- For loga- 
rithmic barriers as in Eq. JSJ, the corresponding growth 
equation is 



dR D 



dt 



R 2 



(1 + R)' e/T . 



(13) 
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The short-time and long-time solutions of Eq. I|13[) are 
obtained as follows: 



R(t) ~ (3£> *) 1/3 , 



where 



e(T,e) = 



0(T,e) 



£> * 



i > to, 



3 + e/T' 

The crossover form of Eq. I|14|) is Eq. © with 



(14) 



(15) 



and 



to = 



h(x) 



(30)0/(1-30) ' 

1 1 

^(3038)1/(1-38): 



(16) 



1/3 , X«l, 

e , a:»l. 



(17) 



Notice that the asymptotic exponent differs from that 
for the nonconserved case when the energy-barriers are 
logarithmic. This should be contrasted with the HH sce- 
nario, where the asymptotic growth law is the same for 
the nonconserved and conserved cases [Io|. This is eas- 
ily seen by incorporating the HH barrier-scaling form in 
Eq. C2J. 



III. RANDOM BOND ISING MODEL: 
CONSERVED KINETICS 

A. Modeling and Numerical Details 

The Hamiltonian for the RBIM is as follows: 

H Y^.J,,S,S r Si = ±l. (18) 

For a binary (AB) mixture, the spins Si label whether 
a lattice site i is occupied by an A-atom (say, Si = +1) 
or a B-atom (Si = —1). We consider the case where the 
spins are placed on an L 2 square lattice with periodic 
boundary conditions. We introduce quenched disorder 
in the exchange coupling , corresponding to immobile 
impurities in a binary mixture. The Jij 's have a uniform 
distribution on the interval [1 — e/2, 1 + e/2], where e 
quantifies the amount of disorder. The limit e — 2 corre- 
sponds to maximum disorder, and e = corresponds to 
the pure case. (We confine ourselves to the case where the 
exchange couplings are always ferromagnetic, Jij > 0.) 
The subscript (ij) in Eq. (|18fl denotes a sum over nearest- 
neighbor pairs only. For the pure case, TP ure ~ 2.269 
for a d = 2 square lattice. Since the average coupling 
strength is (Jij) — 1, as in the pure case, the critical 



temperature remains almost unaltered, T c s [2.0, 2.269] 
|25|. Assigning random initial orientations to each spin, 
we rapidly quench the system to T < T c . The initial 
condition corresponds to a critical quench, with 50 % A 
(up) and 50 % B (down). 

The Ising model has no intrinsic dynamics as the com- 
mutator of the spin variables and the Hamiltonian is iden- 
tically zero. Therefore, we introduce stochastic dynamics 
by placing the system in contact with a heat bath. The 
resultant dynamical model is referred to as a kinetic Ising 
model. The appropriate stochastic kinetics for a binary 
mixture is Kawasaki spin-exchange or conserved kinetics, 
where a randomly-selected spin Si is exchanged with a 
randomly-chosen neighbor, S% <-» Sj. The spin exchange 
is accepted with probability 



W = 



exp(-pAH) for AH > 0, 
1 for AH < 0, 



(19) 



where AH is the change in energy resulting from the spin 
exchange: 

AH = (Si - Sj) ( J > l > Sl > ~ £ Sl ] ' ( 2 °) 



K L i7 tj 



Lj=£i 



In Eq. I|20|l , Li refers to the nearest-neighbors of lattice 
site i. A single Monte Carlo step (MCS) corresponds to 
attempted updates of L 2 spins. A naive implementation 
of the Kawasaki model is numerically demanding, and it 
has proven notoriously difficult to access the asymptotic 
LS growth regime in the pure case [2(| |2tJ • A number 
of accelerated algorithms have been proposed in the lit- 
erature |2S| - we employ the so-called continuous-time 
algorithm. In this approach, a list of oppositely-oriented 
spins is prepared from the lattice configuration. Then, a 
pair is selected randomly from the list, and is exchanged 
according to Eq. (|19fl . In each trial, time is advanced by 
At = l/n t , where nt is the total number of anti-aligned 
spin pairs at time t. After each exchange, the list is up- 
dated. This algorithm works particularly efficiently at 
low temperatures, where bulk domains are strongly en- 
riched in one component. 

The segregating system is usually characterized by 
studying the time-dependence of the correlation function: 



1 N 

C ^t) = nH [(Si(t)S i+f (t)) - (Si(t))(S i+r (t))l v , 
s=i 

. (21) 

which measures the overlap of the spin configuration at 
distance r. Here, [...] av indicates an average over differ- 
ent realizations of the bond disorder, and (...) denotes 
a thermal average, i.e., an average over different initial 
configurations and realizations of the thermal noise. Typ- 
ically, the growth process is isotropic and characterized 
by a unique length scale R(t). In that case, the correla- 
tion function has a dynamical-scaling form p9l | : 



C(r,t)=g 



(22) 
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where g(x) is the scaling function. 

The characteristic size R(t) is defined from the corre- 
lation function as the distance over which it decays to 
(say) zero or half its maximum value. There are a num- 
ber of different definitions of the length scale, but these 
are all equivalent in the scaling regime. Subsequently, we 
will present results for the correlation function and the 
domain growth law. 



B. Numerical Results 

In Fig. ^ we show evolution pictures for the conserved 
RBIM after a critical quench from T = oo to T = 1.0. 
We show snapshots at t = 10 7 MCS for e = (pure 
case), and e = 1,2. The domains have been identified by 



e = 



£= 1 



e = 2 




FIG. 1: Domain growth in the RBIM with Kawasaki kinetics. 
We show evolution pictures a,tt = 10 7 MCS for a 256 2 lattice, 
after a quench from T = oo to T = 1.0. The mixture has 
a critical composition with 50 % A (Si = +1, marked in 
black) and 50 % B (Si = — 1, unmarked). The snapshots 
correspond to different disorder amplitudes: e = (pure case), 
and e = 1,2. 

calculating the time-average for each spin: 



a !>(*)• 



(23) 



within a suitable time- window A = tf — tj. It is clear 
from the snapshots that the evolution is slower for higher 
amplitudes of disorder. This will be quantified via the 
corresponding domain growth laws. 



Next, we consider the scaled correlation-function data 
\C(r,t) vs. r/R] for the morphologies in Fig.^ Our sta- 
tistical data for the RBIM is obtained on d = 2 lattices 
of size 512 2 (with T = 1.0 and e being varied), and 256 2 
(with e = 2 and T being varied). In order to improve 
the statistics, we averaged within a finite time-window 
around each data point. Further, the data was obtained 
as an average over 32 independent initial conditions for 
both the spin and disorder configurations. The length 
scale R is defined as the first zero-crossing of the corre- 
lation function. We have confirmed that C (r, t) exhibits 
dynamical scaling [as in Eq. (|22|l ] for different disorder 
amplitudes and quenches to different values of T. 



O 




r/R 



FIG. 2: Scaling plot of the correlation function for the evo- 
lution depicted in Fig. □ We plot C(r,t) vs. r/R at t = 10 7 
MCS for disorder amplitudes e = (pure case), and e = 1, 2. 
The length scale is defined as the first zero-crossing of C (r, t) . 

In Fig. we show that the scaling function is inde- 
pendent of the disorder amplitude. This has also been 
demonstrated in earlier studies of phase separation in 
disordered systems |ToL Il3|. In physical terms, the uni- 
versality of the scaling function means that the morpholo- 
gies are equivalent, regardless of the disorder amplitude. 
(This was already suggested by the snapshots in Fig.Q]) 
The typical transverse displacement of interfaces due to 
disorder roughening is L^ 2- ^, where £ is the roughen- 
ing exponent |5(. At late times, one has L > L^ 2 -C) 
because £ < 1 above the lower critical dimension. (If 
C > 1, disorder-induced roughening would destroy long- 
range order in the system.) Thus, in the asymptotic 
regime, the roughness is irrelevant compared to the do- 
main size. Therefore, the evolution morphologies and 
their statistical properties should be independent of dis- 
order at late times. 

Next, let us investigate the time-dependence of the do- 
main size. First, we study R(t) vs. t for quenches to 
different temperatures. 

In Fig. |3 we undertake a direct test of the HH growth 
law in Eq. © by plotting _R 1/4 vs. baf for e = 2 and dif- 
ferent T- values. Recall that ip = 1/4 in d = 2 according 
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FIG. 3: Plot of i? 1 / 4 vs. t (on a log-linear scale) for e = 2 and 
different quench temperatures: T — 0.7, 0.8, 0.9, 1.0, 1.2. 



to the HH argument, and the corresponding plot in Fig.|3 
should be linear in the asymptotic regime. However, the 
plot exhibits continuous curvature and is not consistent 
with the HH growth law. We have also attempted to fit 
the data to the functional form lnt = aR x + b. In general, 
this function does not give a reasonable fit to the data. 
Even for these poor fits, the exponent x is strongly depen- 
dent on the temperature, at variance with the prediction 
of a universal growth law. A similar observation has been 
made in the experiments of Ikeda et al. [18( , though these 
were performed on random magnets, rather than disor- 
dered mixtures. As a matter of fact, Ikeda et al. and 
Likodimos et al. 

El Hi 

have argued that their experi- 
mental data for domain growth in disordered systems is 
described by a power-law with a temperature-dependent 
exponent rather than the HH growth law. We have made 
a similar observation in our MC studies of the noncon- 
served RBIM |l4j. Let us examine the length-scale data 
for the conserved RBIM from this perspective. 

In Fig. 0Ja), we plot R vs. t from Fig. [3] on a log- log 
scale. This plot does not show an extended linear regime 
on the time-scale of our simulation. However, it is known 
that there is an extended pre-asymptotic grow th regime 
in the conserved case without disorder |24l l26l |27| , which 
complicates the observation of the LS growth regime in 
MC simulations. Further, the slight upward curvature 
in the log-log plot suggests that the growth law cannot 
be slower than a power law, at variance with the HH re- 
sult. In the pure case, Huse [24[ has suggested that the 
asymptotic exponent may be obtained by extrapolating 
the graph of the effective exponent 9 c g — d(\n R) / d(\n t) 
vs. Rr 1 . We apply a similar technique to the disordered 
case, and query whether the resultant exponents are con- 
sistent with the scenario in Sec. II. B [cf. Eq. I|15|) ]. In 
Fig. Etb), we plot 9 cS vs. R^ 1 for the data in Fig.^Ja). 
The plots in Fig. Efb) can be smoothly extrapolated to 
Rr 1 = (R = oo ) to determine 6 = 9 e g(oo), which de- 
pends on T. 
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0.4 
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(b) 




0.1 



0.2 
1/R 



0.3 



0.4 



FIG. 4: (a) Plot of R vs. t (on a log-log scale) for the length- 
scale data shown in Fig. [3] (b) Plot of 8 e s = d(ln R) / d(ln t) 
vs. _R _1 for the data in (a). 



Next, we consider R vs. t at fixed temperature as 
the disorder amplitude is varied. Again, we find that our 
data is not consistent with either the HH scenario or even 
logarithmic growth. In Fig.OJa), we plot R vs. t on a log- 
log scale for different e-values. The corresponding plots 
of 9 e s vs. are shown in Fig. |S£b). In this case, the 
asymptotic exponent depends on the disorder amplitude. 
Notice that we have also shown data for the pure case 
(e = 0) in Fig. [SJa). We do not see an extended linear 
regime even in this case. However, the corresponding 
plot of 8 c s vs. R^ 1 in Fig. EJb) extrapolates to the well- 
known LS value, 9 ~ 0.33. 

In Sec. II. B, we have seen that a logarithmic barrier- 
scaling results in power-law growth with varying expo- 
nents. We would like to test whether the asymptotic 
exponents are consistent with the result in Eq. (|15f) . In 
Figs. EJa) and (b), we plot 9~ 1 vs. T _1 and e, respec- 
tively. The resultant linear plots strongly support the 
logarithmic barrier-scaling scenario. 
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0.1 1 1 1 1 1 

0.1 0.2 0.3 0.4 

1/R 

FIG. 5: (a) Plot of R vs. t (on a log-log scale) for T = 1.0 
and different disorder amplitudes: e = (pure case), and 
e = 0.5, 1.0, 1.5, 2.0. (b) Plot of 6 cS = d(]nR)/d(]nt) vs. R' 1 
for the data in (a). 



3®^- — ■ ■ ■ 

0.5 1 1.5 2 

e 

FIG. 6: (a) Exponent 1/6 vs. l/T for the data in Fig. The 
solid line denotes the best linear fit to the data, (b) Exponent 
1/6 vs. e for the data in Fig. EI 



IV. DILUTE ISING MODEL 

A. Modeling and Numerical Details 

Next, we turn our attention to the DIM, where bond 
disorder is introduced by diluting the spins on the lattice. 
The corresponding Hamiltonian is 

H = -J^e i E j S i S j , Si=±l, (24) 

(ij) 

with J > 0. In Eq. Ifity . the e^'s are quenched, uncorre- 
cted random variables with the probability distribution: 

P(e)=pS e>1 + (l-p)S et0 . (25) 

For a ferromagnet, = implies that the magnetic atom 
at i is replaced by a non-magnetic impurity. In the con- 
text of an AB mixture, Ei — corresponds to an immo- 
bile (non-interacting) impurity at site i. Thus, there is 
no exchange interaction between the atom at site i and 
its nearest neighbors. The distinguishing feature of the 



DIM (in contrast to the RBIM discussed in Sec. Ill) is 
the existence of a percolation threshold p = p c [3(j • For 
p = 1, the system is pure and shows ferromagnetic or- 
der at T < T c (p = 1). The critical temperature T c (p) 
diminishes as p is decreased and becomes at p = p c . 
(For a d — 2 square lattice, p c ~ 0.593.) For p < p c , 
there are no infinite clusters of magnetic atoms which 
span the system, i.e., there is no long-range order. For 
weak disorder (p ~ 1), the kinetic DIM is analogous to 
the kinetic RBIM. However, for smaller values of p, con- 
nectivity effects become important and may change the 
nature of domain growth. We are particularly interested 
in the ordering dynamics of the DIM for p ~ p c . 

In this section, we focus on two systems: 

(a) The DIM with nonconserved (Glauber) kinetics, 
which models the ordering dynamics of a dilute ferro- 
magnet. In an MC simulation of Glauber kinetics, a 
randomly-chosen spin Si is flipped to — Si and the system 
is evolved according to the prescription in Eq. (|19|) . 

(b) The DIM with conserved (Kawasaki) kinetics, which 
models the segregation kinetics of a dilute binary mix- 
ture. In this case, we use the continuous-time algorithm 
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described in Sec III. A. 

The initial conditions for our MC simulations are pre- 
pared as follows. We dilute the sites of an 1? lattice 
with probability 1 — p. (These sites remain fixed during 
the evolution.) Then, up and down spins are randomly 
distributed on the remaining sites with a zero net mag- 
netization, mimicking the high-temperature disordered 
configuration before the quench. 



B. Nonconserved Kinetics 

In Fig. 13 we show evolution snapshots at t = 10 6 MCS 
for T = 0.5 and p = 0.9, 0.8, 0.7, 0.593 (p c ). Notice that 



mm, 




FIG. 7: Domain growth in the DIM with Glauber kinetics. 
We show evolution pictures at t = 10 6 MCS for a 256 2 corner 
of a 512 2 lattice, after a quench from T = ootoT = 0.5. The 
snapshots correspond to different site occupation probabili- 
ties: p = 0.9,0.8,0.7,0.593 (p c ). The up and down spins are 
marked black and white, respectively. The missing spins are 
marked grey. 

T c (p = 0.7) ~ 1.04 for the d = 2 DIM H3, so that 
T = 0.5 lies below the critical temperature for all the 
values of p other than p = p c , where T c (p c ) = 0. (Unfor- 
tunately, it is difficult to do MC simulations at T = 0, 
as the system is rapidly trapped in a metastable state.) 
As expected, the domain size at a fixed time diminishes 
with increase in disorder. In the case of evolution on 
the backbone of a percolating cluster, the morphology 
consists of a network of islands (compact well-connected 
regions) linked by just a single bond. As time progresses, 
these islands become fully magnetized, but they cannot 



influence the evolution of their neighbors. For this rea- 
son, domain growth becomes very slow at p = p c . Fur- 
ther, as T > T c (p c ), domain growth is arrested when 
the length scale saturates at the equilibrium correlation 
length £ eq (T) -> oo as T -► 0. 




r/R 



FIG. 8: Scaling plot of the correlation function for the evo- 
lution depicted in Fig. □ We plot C(r,t) vs. r/R at t = 10 6 
MCS for occupation probability p — 0.8, 0.7, p c . We also show 
data for the pure case (p = 1) at t — 10 3 MCS, obtained for 
a 1024 2 system. 

Let us next focus on the properties of these evolution 
morphologies. Our statistical data for the nonconserved 
DIM is obtained using 512 2 systems, by averaging over 
50 independent initial conditions and disorder configu- 
rations. We have confirmed that the evolution of the 
nonconserved DIM shows dynamical scaling. In Fig. [SI 
we demonstrate the disorder-independence of the scaled 
correlation function. Here, we plot C(r, i) vs. r/R at 
t = 10 6 MCS for p = 0.8, 0.7, p c , and compare it with 
the corresponding data for the pure case (p = 1). In this 
case, the domain size is defined as the r- value where the 
correlation function decays to half its maximum value. 
Notice that the scaling function for p = p c is analogous 
to that for higher values of p, and there are no distinctive 
signatures of the percolation cluster. 

Next, consider the time-dependence of the length scale. 
We first study the case with p = 0.8 and varying T- values. 
In Fig. El we attempt to fit our length-scale data to the 
HH crossover function in Eqs. |JS}-|J5J. We record the 
following points of disagreement with the HH scaling be- 
havior: 

(a) The short-time behavior is not described well by 
Eq. where h\{x) ~ x 1 / 2 . 

(b) The asymptotic behavior in Eq. JSJ [denoted by the 
curve h,2 (x) in Fig. [5] does not fit the scaling curve well 
even for the largest times. 

(c) The temperature-dependence of the crossover length 
Rq(T) and the crossover time to(T) is stronger than a 
power law (see inset), which is incompatible with Eqs. @ 
and J3J|. The parameter ao in Eq. is proportional to 
the surface tension, and is expected to decrease with in- 
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FIG. 9: Scaling plot to test the crossover function in Eqs. JjJJ- 
©. For each temperature T, the values for Ro(T) and to CO 
have been chosen to obtain a smooth scaling curve h{x). The 
functions h\(x) oc a; 1 ' 2 and Ii2(x) oc (lnx) 4 represent the 
expected asymptotic behavior for i < 1 and respec- 
tively. The inset shows the temperature-dependence of the fit 
values Ro(T) and to(T), and their expected T-dependence, 
which is T 4 and T 8 , respectively. 



creasing temperature. Therefore, t can be expected to 
increase faster than T 8 , but its T-dependence turns out 
to be much too strong: note that to in the inset of Fig. [5] 
varies over 20 decades when T varies over only half a 
decade from 0.4 to 0.8. We do not see why the surface 
tension should have such a strong T-dependence. 
Based on observations (a)-(c), we believe that the data 
in Fig. is inconsistent with the HH growth law. In 
Fig. HUT a). we plot R vs. t on a log- log scale for p = 0.8 
and T = 0.8,0.7,0.6,0.5,0.4 < T c (p = 0.8) ~ 1.5. 
The corresponding plots of 9 c tf vs. R~ x are shown in 
Fig. ITUT bh These show an extended flat regime, making 
it relatively simple to estimate the exponent. As in the 
case of the RBIM, our data is consistent with power-law 
growth with a variable exponent. In Fig. HOf ch we plot 
#(T,p) _1 vs. T _1 - the linear behavior is consistent with 
Eq. ©. (See Ref. for similar results for the noncon- 
served RBIM.) 

Finally, in Fig. Illf ah we plot R vs. t at percolation 
(p = p c ) and T = 0.7,0.6,0.5,0.4 > T c (p c ) = 0. Recall 
that the domain scale saturates to £ cq (T) in this case, 
with an earlier crossover for higher T. On the time-scale 
of our simulation, the data for T = 0.7 has saturated, 
and that for T = 0.6 is beginning to bend over. This is 
reflected in Fig. Illf b). which shows 9 e g vs. R . The 
exponent 9 is estimated from the flat portion of these 
curves, and we plot 6* _1 vs. T _1 in Fig. Illf ch 
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FIG. 10: (a) Plot of R vs. t (on a log-log scale) for the 
nonconserved DIM with p — 0.8 and temperatures T — 
0.4,0.5,0.6,0.7,0.8. (b) Plot of cft • = d(ln R)/d(\n t) vs. iT 1 
for the data in (a), (c) Plot of vs. T _1 for the data in 
(a). 



C. Conserved Kinetics 

We have performed a similar study of the DIM with 
Kawasaki kinetics. In this case, the time-scale of growth 



is considerably slower than for the nonconserved case. 
The typical evolution morphologies at t = 10 7 MCS (after 
a critical quench from T = ootoT = 0.5) are shown in 
Fig. E| As in the earlier cases, we will show results 
for the correlation function and the growth law. The 
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FIG. 11: Analogous to Fig.^but for p = p c 



statistical data shown here was obtained on a 256 2 lattice 
as an average over 32 independent configurations. 

In Fig.HBI we plot C(r,t) vs. r/R at t = 10 7 MCS for 
the pure case, and different values of the dilution. [The 
length scale is obtained from the first zero of C(r, t).] 
Again, the scaling function is approximately independent 
of the amount of dilution. Next, we focus on the time- 
dependence of the length scale. In Fig. I14f a). we plot 
R vs. t for p — 0.8 and various values of T. Again, 




p = 0.7 



'if*. 




FIG. 12: Domain growth in the DIM with Kawasaki kinetics. 
We show evolution pictures at t = 10 7 MCS for a 128 2 corner 
of a 256 2 lattice, after a quench from T — oo to T = 0.5. The 
snapshots correspond to different site occupation probabili- 
ties: p — 0.9, 0.8, 0.7, 0.593 (p c ). The color coding is the same 
as in Fig. Q 




FIG. 13: Scaling plot of the correlation function for the 
evolution depicted in Fig. 1121 We plot C(r,t) vs. r/R at 
t = 10 7 MCS for occupation probability p — 1 (pure case) 
and 0.8, 0.7, p c . 



we estimate the asymptotic exponent from plots of 9 c g 
vs. R^ 1 (not shown here). In Fig. I14f b). we plot the 
corresponding 6" 1 vs. T _1 . 

FigurelTSlis analogous to Fig.^1 but for p = p c . As the 
growth is much slower than the nonconserved case, we do 
not see a crossover to saturation for p = p c on the time- 
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FIG. 14: (a) Plot of R vs. t (on a log-log scale) for the 
conserved DIM with p = 0.8 and T = 0.4, 0.5, 0.6, 0.7, 0.8. (b) 
Plot of 6» _1 vs. T" 1 for the data in (a). 



scale of our simulations. Once again, the exponents are 
consistent with the logarithmic barrier-scaling scenario. 



V. SUMMARY AND DISCUSSION 

Let us conclude this paper with a summary and dis- 
cussion of the results presented here and in our earlier 
letter 14]. We have undertaken comprehensive Monte 
Carlo (MC) simulations of domain growth in Ising sys- 
tems with quenched disorder. These studies are based on 
kinetic Ising models with either nonconserved (Glauber) 
spin-flip kinetics or conserved (Kawasaki) spin-exchange 
kinetics. The nonconserved case models ordering dynam- 
ics in random magnets, and the conserved case models 
segregation kinetics in disordered binary mixtures. We 
have studied domain growth for two classes of disordered 
systems: 

(a) The random-bond Ising model (RBIM), where the ex- 
change interaction has a uniform distribution on the in- 
terval [1 — e/2, 1 + e/2], e < 2. In this case, the critical 
temperature T c (e) remains approximately unchanged. 



FIG. 15: Analogous to Fig. 1141 but for p = p c . 



(b) The dilute Ising model (DIM), where the exchange 
interaction is randomized by the dilution of magnetic 
atoms with non-magnetic impurities. In this case, the 
critical temperature T c (p) ranges from T c (p = 1) ~ 2.269 
(in d = 2) to T c (p = Pc )=0 ( Pc ~ 0.593 in d = 2). 
Both classes of disorder are of considerable experimental 
relevance. 

The general framework for understanding coarsening 
in disordered systems is as follows. At early times, the 
domain sizes are small and domain growth is unaffected 
by disorder. At late times, the domain boundaries are 
trapped by disorder sites, and asymptotic growth pro- 
ceeds via thermally-activated hopping over disorder bar- 
riers. Clearly, the asymptotic growth law depends criti- 
cally on the length-scale dependence of the disorder bar- 
rier Eb- In this context, an important study is due to 
Huse and Henley (HH) 0. In the HH scenario, the disor- 
der barriers have a power-law dependence on the domain 
size, Eb ~ ■ These result in a logarithmic domain 
growth law in the asymptotic regime. We find that our 
MC results are not in agreement with the HH scenario. 
Rather, our results are consistent with power-law growth 
with an exponent 9 which depends on the temperature 
T and the disorder amplitude e. This is in agreement 
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with a number of experiments 0, 0, |2(|, and early 
simulations of droplet shrinking in disordered systems 
by Oh and Choi Q- This scenario arises naturally in 
the context of logarithmic energy barriers, and the corre- 
sponding functional dependence of 0(T, e) is in excellent 
agreement with our numerical results. 

Our results provide a framework for the analysis of 
experiments and simulations on domain growth in dis- 



ordered magnets and binary mixtures. We hope that 
our study will motivate fresh experimental studies of this 
important problem. In particular, there is a paucity of 
experimental results on phase separation in disordered 
mixtures. 
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